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ABSTRACT 


A  compact  heat  exchanger  testing  facility  was  constructed  at  the 
United  States  Naval  Postgraduate  School.  The  testing  facility  is  used 
to  collect  heat  transfer  and  friction  factor  data  by  the  steady  state 
steam  to  air  technique.  Heat  exchangers  with  frontal  dimensions  of  up 
to  12  inches  by  12  inches  can  be  tested.  The  maximum  pressure  drop  which 
can  be  produced  across  the  air  side  of  a  heat  exchanger  is  2.5  psi.  The 
testing  facility  was  designed  to  obtain  data  for  exchanger  air  side 
Reynolds  numbers  from  500  tc  10,000.  A  description  of  the  apparatus,  the 
experimental  method  and  equations  and  an  evaluation  are  given. 
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NOMENCLATURE 


English 
Aa 
Ac 

Afa 
Afr 

Af  s 

As 

^va 

^ws 
r i 


Letter  Symbols 

2 

-  Air  side  total  heat  transfer  area,  ft 

-  Air  side  free  flov  area,  ft 

-  Air  side  fin  area,  ft^ 

2 

-  Air  side  frontal  area,  ft 

-  Steam  side  fin  area,  ft2 

-  Steam  side  total  heat  transfer  area,  ft2 

-  Air  side  wall  area,  ft2 

-  Steam  side  wall  area,  ft2 

-  Coefficient  of  discharge  for  a  metering  orifice 

Cc  -  Jet  contraction  ratio 

Cg  -  Coefficient  of  discharge  of  a  fluid  orifice  for  a  given  value 
of  3 

-  Specific  heat  of  air  at  constant  pressure,  BTU/(lbm  °F) 

d  -  Diameter  of  fluid  metering  orifice,  in. 

D  -  Air  orifice  duct  diameter,  in. 

D<j  -  Hydraulic  diameter  of  duct  immediately  downstream  of  test 
core,  in. 

ERROR  -  Error  in  heat  balance  between  steam  side  and  air  side,  percent 
f  -  fanning  friction  factor  in  test  core 

f<j  -  Friction  factor  in  duct  immediately  downstream  of  test  core 

fm  -  Friction  factor  of  a  passage  in  test  core 

F  -  Velocity  of  approach  factor  for  the  fluid  metering  orifice 

Fa  -  Thermal  expansion  factor  for  the  fluid  metering  orifice 

gc  -  Proportionality  factor  in  Newton's  Second  Law,  gc  ■  32.174 
lbm  ft/(lbf  sec2) 

6  -  Test  core  air  mass  velocity,  !bsm/ (hr  ft2) 
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Unit  conductance  for  thermal  convection  heat  transfer 
(film  coefficient),  BTU/(hr  ft2  °F) 

Enthalpy  of  the  condensate  leaving  the  test  core,  BTU/lba 

Enthalpy  of  steam,  BTU/lbm 

Humidity  ratio  for  air,  Ibjn  water/lbm  dry  air 

Unit  thermal  conductivity,  BTU/ (hr  fi:2  °P/ft) 

Contraction  coefficient 
Velocity  distribution  coefficient 
Expansion  coefficient 

Fin  length  equal  to^one  -half? of- total  fin  height 

Length  of  duct  from  downstream  face  of  test  core  to  downstream 
pressure  tap,  in. 

Wall  thickness  between  steam  and  air  side  of  test  core,  in. 
Test  core  flow  length,  ft. 

VTh7k7 

Mass  rate  of  air  flow,  Ib^/hr. 

Pressure,  in.  H2O,  in.  Hg,  lbf/in2,  lbf/ft2 

Barometric  pressure,  in.  Hg 

Heat  transfer  rate,  BTU/hr 

Gas  constant  for  air,  53.3  ft  lbf/(lbm  °R) 

Hydraulic  radius,  (A^/A^),  ft,  (4  rj,-hydraulic  diameter) 

Mass  rate  of  steam  flow,  lb  m/hr 

Temperature,  °F 

Absolute  temperature,  °R,  °K 

Unit  overall  thermal  conductance,  BTU/(hr  °F  ft2) 

Velocity,  ft/sec 
Specific  volume,  ft2/lbm 

Mass  rate  of  condensate  from  test  core,  loi^/hr 


X 


Xm 

Y 


-  Ratio  of  pressure  upstream  of  the  fluid  metering  orifice  to 
pressure  differential  across  the  orifice 

-  Humidity  correction  to  the  specific  heat  of  air 

-  Humidity  correction  to  the  density  of  air 

-  Net  expansion  factor  for  a  square-edged  metering  orifice 


Greek  Letter  Symbols 

5  -  Ratio  of  fluid  metering  orifice  diameter  to  duct  diameter 

A  -  Denotes  difference 

6  -  Fin  thickness,  in. 

n  -  Fin  temperature  effectiveness 

n0  -  Total  surface  temperature  effectiveness  for  air  side 

i13  -  Total  surface  temperature  effectiveness  for  steam  side 

<T  -  Rat.'.o  of  free-flow  to  frontal  area,  A’/AfR 

/t  -  Dynamic  viscosity,  lbm/ (hr  ft) 

p  -  Density,  lb  -fa/ft3 

"  Average  density  of  air  in  test  core  for  isotherma]  test, 

lb  B/ft3 


% 


-  Mean  density  of  air  in  test  core  for  hot  core  test,  lb.*m/ft“ 


Dimensionless  Groupings 

%e  ~  Reynolds  number,  (Ar^G^) 

Ngt  -  Stanton  number,  (h/GCp) 

NPr  -  Prandtl  number,  ^/UCp/k) 

j  -  Colburn-j  (NSt  Npr2/3) 

NTU  -  Number  of  heat  transfer  units  for  an  exchanger  (Aa  U/fi  Cj) 
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Subscripts 

1  -  Air  upstream  of  core 

2  -  Air  downstream  of  core 

o  -  Air  metering  orifice 

a  -  Air  side  of  test  core 

c  -  Core 

f  -  Fin 

f  -  Air  film  in  test  core 

s  -  Steam  side  of  test  core 

so  -  Steam  metering  orifice 

si  -  Steam  in  top  header 

s2  -  Steam  in  bottom  header 

w  -  Wall  between  air  and  steam  sides  of  test  core 

s  -  Square  fin 

t  -  Triangular  fin 

c  -  Circular  tube 

Force  and  Mass  Units 

lbm  -  Denotes  pounds  mass  in  distinction  to 

lbf  -  Denoting  pounds  force 
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1.  Introduction 


Steady  state  tests  transferring  heat  from  condensing  steam  to 
flowing  air  are  used  to  determine  experimentally  the  heat  transfer 
characteristics  for  compact  heat  exchangers.  The  friction  factor  is 
determined  from  the  air  pressure  drop  across  the  heat  exchanger.  Steam 
and  air  are  used  as  the  working  fluids  in  the  tests  because  they  are 
readily  available  and  produce  effects  which  are  easily  analyzed.  The 
heat  transfer  characteristics  are  given  by  Colburn-j  heat  transfer 
modulus  evaluated  from  the  steady  state  steam  to  air  test  (referred  to 
here-in  as  hot  core  test).  The  friction  factor  is  evaluated  from  both 
the  hot  core  test  and  the  isothermal  test  (referred  to  as  cold  core  test) 
The  Colburn-j  heat  transfer  modulus  and  friction  factor  are  normally 
presented  graphically  as  functions  of  the  heat  exchanger  Reynolds  number. 
Development  of  these  three  parameters  is  given  in  the  theory  section. 

In  addition  to  determining  the  Colburn-j  and  friction  factor.,  a 
heat  balance  is  performed  on  the  exchanger  to  check  the  consistency  of 
the  experimental  data. 

The  testing  facility  is  similar  in  design  to  one  built  at  Stanford 
University  in  1947.  The  Stanford  facility  was  built  under  a  United 
States  Navy  contract  to  evaluate  finned  surfaces  for  use  in  compact  hdat 
exchangers  and  is  described  by  Kays  in  reference  [3]. 

Compact  heat  exchangers  which  are  suitable  for  steady  state  steam 
to  air  testing  have  plate-fin  geometry  as.  shown,  below. 
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The  steady  state  steam  to  air  testing  facility  was  designed  and 
built  to  accommodate  a  heat  exchanger  with  frontal  dimensions  of  up  to 
12  inches  by  12  inches.  The  facility  is  capable  of  producing  air  side 
Reynolds  numbers  from  500  to  10,000.  An  air  side  pressure  drop  of  2.5 
psi  can  be  obtained. 

The  testing  facility  is  described  and  the  experimental  method  and 
equations  for  data  reduction  are  presented.  The  data  reduction  equations 
are  conveniently  arranged  for  programing  on  a  digital  computer.  A  test 
of  a  heat  exchanger  of  known  characteristics  is  presented  as  an  evalua¬ 
tion  of  the  testing  facility.  The  operating  procedures  and  the  computer 
program  for  data  reduction  with  accompanying  instructions  are  given  in 
the  appendix. 
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2.  Theory 


Heat  Transfer  Characteristics.  The  heat  transfer  characteristics 
of  a  compact  heat  exchanger  are  best  represented  by  dimensionless  group¬ 
ings.  This  generalizes  the  results  so  that  they  may  be  used  for  working 
fluids  ether  than  steam  and  air.  The  succeeding  development  is  that 
commonly  used  in  the  field;  see  Kays  and  London  [4].  The  Colbura-j  heat 
transfer  modulus  is  selected  for  this  purpose  as  it  combines  heat  transfer 
characteristics  and  fluid  properties.  The  Colburn-j  modulus  is  defined  as: 


j  5  ^St  NPr 


G  Cr 


Jisz, 
.  k  ) 


y/3 


where  the  .unit  conductance  for  convective  heat  transfer,  b,  must  be  deter¬ 
mined  experimentally.  The  use  of  a  steam  to  air  system  is  one  of  the 
most  direct  methods  of  experimentally  determining  h. 

The  unit  conductance  for  convective  heat  transfer,  h,  for  the  air 
side,  air  being  the  fluid  of  interest,  is  determined  from  energy  balance 
considerations  within  the  exchanger.  Low  pressure  slight  .y  superheated 
steam,  about  6  psig  with  5  degrees  superheat,  enters  the  top  of  the  heat 
exchanger  steam  side.  The  steam  reaches  its  saturation  temperature  in  a 
very  short  distance  after  entering  the  exchanger.  Constant  temperature 
is  maintained  as  the  steam  condenses  for  the  remainder  of  travel.  The 
air  is  heated  continuously  as  it  passes  through  the  exchanger.  The 
sketch  below  illustrates  the  temperature  variations  in  the  working  fluids 
as  a  function  of  the  heat  transfer  area. 


An  energy  balance  for  the  differential  area  dA  yields  the  equation 
for  the  heat  transfer  rate  as: 

dq  =*  1  Cp  dTair  *  U  d  A  AT 

(heat  transfer  (heat  transfer  rate  (heat  transfer  rate 
rate)  to  the  air)  from  steam  to  air) 

Rearrangement  gives: 

-  dTair _  «  _  U  d  A 

(Tsteam-Tair)  m  Cp 

For  small  over-all  changes  in  fluid  properties  and  a  relatively  uniform  flow 
distribution  in  the  heat  exchanger,  the  over-all  unit  conductance,  U, 
may  be  assumed  to  be  constant.  Integration  over  the  total  area  then 
yields: 

Tsteam  -  Tair  out 
Tsteam  -  Tair  in 

The  dimensionless  grouping  on  the  right  side  of  the  equation  is  defined 
as  the  Number  of  heat  Transfer  Units,  NTU.  Thus  the  equation  may  be 
written  as: 


in  I  Tsteam  -  Tair  out 
Tsteam  -  Tair  in 


or 


(- 


Tsteam  -  Tair  out 
Tsteam  -  Tair  in 


-) 


-  NTU 


-NTU 


The  over-all  unit  conductance  based  upon  the  total  air  side  heat 
transfer  area  must  include  a  factor  to  account  for  the  temperature 
variation  in  the  fin  areas.  The  over-all  temperature  effectiveness,  n, 

may  be  written  as:  ^ 

n  =  1  -  (1  -Hf) 

A 
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where  Af  is  the  fin  area  andT^f  is  the  fin  temperature  effectiveness. 
The  effectiveness  for  a  fin  of  constant  conduction  cross  section  as 
given  by  Kays  and  London  [A]  is: 


tanh  ml 
nf  "  ml 


where 


m 


for  thin  sheet  fins 


1  ■  one  half  fin  height  for  wall  to  wall  fins 


Summing  the  resistances  to  heat  flow  then  yields: 


1  _ Ah  i  Aa  x  _ 1 

U  “  nbh 


(over-all 

resistance) 


(steam  side 
resistance) 


(wall 

resistance) 


(air  side 
resistance) 


The  convective  coefficient  for  the  air  side,  h,  may  be  solved  for  giving: 


^  U--_Ad(-n,Ashs  +  A ft)' 


(1) 


A  value  of  hg  equal  to  2000  BTU/ (hr  ft^  °F)  is  assumed.  A  fifty 
percent  error  in  this  assumption  results  in  approximately  two  percent 
error  in  the  value  of  h.  Accordingly  the  heat  transfe."  coefficient  for 
the  air  side,  h,  may  be  evaluated  by  experimentally  measuring  the  fluid 
temperatures  and  the  air  mass  rate  of  flow. 

Friction  Factor.  The  flow  of  a  fluid  through  an  exchanger  results 


in  a  pressure  differential  caused  by  the  flow  contraction  as  the  entrance, 
the  fluid  friction  in  the  exchanger  and  the  flow  expansion  at  the  exchanger 


exit.  These  pressure  effects  are  graphically  represented  below. 


HEAT 

1  EXCHANGER  2 


To  define  a  meaningful  heat  exchanger  friction  factor,  the  entrance 
and  exit  pressure  effects  must  be  elimated  from  the  total  pressure 
differential.  The  friction  factor  will  then  be  based  solely  upon  the 
mean  wall  shear.  The  problem  of  entrance  and  exit  pressure  effects  has 
been  investigated  both  analytically  and  experimentally  by  Kays. [2]. 

The  effect  of  pressure  loss  in  the  downstream  duct  is  included  in 
the  total  pressure  drop  in  order  that  point  2  may  be  located  downstream 
sufficiently  to  assure  complete  pressure  recovery.  The  total  pressure 
differential  from  point  1  to  point  2  is  then: 

p total  "  +  AP'c  -  ^b  +  APabct 

The  expression  for  the  pressure  effect  due  to  flow  contraction  at 
the  entrance  is  given  by  Kays  f2]  is: 

A pd  _  yL(i  "  <rz)  .  KcvJ 

?«3  2gc  \  2  gc 
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and  for  Che  flow  expansion  ac  the  exit  Is: 

Aph  _  v|_(i  -  <rz) 

Pb  “  2gc 


mL 

2gc 


with  <J~  denoting  the  exchanger  free  flow  to  frontal  area  ratio.  Kc  And 
Ke  are  contraction  and  expansion  coefficients  which  are  functions  of 
exchanger  geometry  and  Reynolds  number.  Noting  that  the  velocity  may  be 
related  to  the  mass  velocity  and  that  the  density  change  from  section  one 
to  a  and  b  to  two  will  be  very  small  compared  to  the  total  density  change, 
the  above  equations  yield: 

.2 


Ap  -  6*-~ 

d'2gc?i 


(i-  or2 


Kc) 


and 


Apb  =  I^  (l-ff-2  -  Ke) 


The  pressure  drop  in  the  exchanger  may  be  found  by  writing  Bernoulli’s 
equation  in  the  differential  form  to  include  the  Fanning  friction  term. 
Using  the  mass  velocity,  6,  the  equation  as  given  by  Kays  [3]  is: 


_dPc  G41 

<?  "  zgt 


,  -  a2  (dx 

+  r2gl?  lrh 


or 


,( i  \ 

w) 

62  L  (_1_\  dx  1 

2gc  r\e/rH  ?2  J 

.  the 

f- 


V  - 

m  \?mi  ~  L 


By  defining  the  mean  specific  volume  in  the  exchanger  as: 

-  ±  \ 

e 

and  integrating  from  p^,  0  to  p\j,  L,  the  exchanger  pressure  differential 
is  given  by: 

GZ 

- *1  ~rnr\~l  , 

•my 


Af>c=  2gcPa 


■i 


with  LMTD  being  the  Log  Mean  Temperature  Difference. 


The  pressure  loss  in  the  downstream  duct  may  be  evaluated  as: 

Abduct  -  4  d  \  v<? 

\  Dd  /  2gc 

where  the  density,  p<j>  the  friction  factor,  f^,  the  length  of  duct 
from  section  b  to  2,  the  hydraulic  diameter,  D,j,  and  the  velocity, 
V<j,  are  all  quantities  associated  with  the  downstream  duct.  The  velo¬ 
city  in  the  downstream  duct  may  be  related  to  quantities  appearing  in 
the  previous  pressure  loss  equations.  If  the  downstream  duct  has  the 
same  dimensions  as  the  exchanger  and  small  density  variations  are 
neglected,  the  continuity  equation  gives 


or 


Thence 


Vd 

Vd 


/ft  )[  I  \('Ac\ 

\Ac/l  ?b)\  Arr/ 


.S-JL. 

?b 


APr 


duct 


=  4 


G2  C~2 

^9c  Cb 


Again  assuming  that  the  density  changes  between  section  1  to  a  and 


b  to  2  are  negligible,  the  various  pressure  drops  are  summed  yielding: 

'z  r(i-<r2+Kc)  2 (Pi  -^2) 

.  e,  + 


AR  t,.  =-§_ 
total  2gc 


rhPm 

Manipulation  of  the  equation  gives: 


■  fL  _(l-(T2-K»)  ,  4flk?i|(r‘ 

+  rhP„  V - 


Aptotal  =  2^ 


fL 

rhe, 


m 


+ 


T2 

Kc-(cr2t-i) 

—?r~ 


1 


and  solving  the  equation  for  the  friction  factor  we  obtain: 


f  = 


£h&n 

L 


ptofal 

G2/2gc 


(kc~  cr2 -I) 

P, 


-  £.+«,+ c=fn-M  ■ 

?2  J 

Therefore  the  friction  factor,  f,  for  a  heat  exchanger  of  given  geometry 
may  be  found  by  experimentally  determining  the  mass  rate  of  air  flow, 

the  entrance  and  exit  air  pressures  and  temperatures  and  the  steam 
temperature. 


(2) 


■Rsyndds  Number*  For  uniformity  of  reporting  test  results,  the 
Reynolds  number  is  defined  in  the  same  manner  as  given  by  Kays  and 
London. [4].  This  definition  embodies  the  hydraulic  diameter  as  the 
characteristic  length  and  the  velocity  based  cn  the  minimum  free  flow 

(\ 

area.  Thus:  1“ 

*Re  -  ^  =  4p£~  =  ^  • 

M  M  ac  /t 

* 

To  evaluate  the  foregoing  characteristics,  the  test  apparatus 
must  be  capable  of  accurately  measuring  the  air  state  upstream  and 
downstream  of  the  exchanger,  the  steam  state  within  the  exchanger  and 
the  mass  rate  of  air  flow. 


3.  Description  of  the  Test  Apparatus 

£aoeraA.  Air  is  induced  through  one  side  of  the  compact  heat 
exchanger  by  a  centrifugal  compressor  and  associated  ducting.  Steam 
is  passed  through  the  other  side  of  the  exchanger  by  the  steam  system. 
The  air  and  steam  systems  are  shown  diagramatically  in  Fig.  4  and  Fig.  5. 
The  steam  system  is  designed  to  give  close  control  of  the  steam  state 
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within  the-  core,  pressure  and  temperature  instrumentation  are  pro¬ 
vided  for  measuring  the  air  properties  upstream  and  downstream  -of  the 
core,  the  mass  rate  of  air  flow,  the  condition  of  the  steam  entering 
and  leaving  the  core  and  the  mass  rate  of  excess  steam  and  condensate 
leaving  the  core. 

The  test  apparatus  is  shown  pictorially.  in  Figs.  1,  2,  and  3. 

Air  System.  The  air  system  is  made  cp  of  16  gage  galvanized 
steel  ducts  joined  with  1/4  inch  steel  flanges.  The  ducting  prior  to 
the  test  section  has  a  maximum  dimension  of  12  inches  square  and  the 
ducting  after  the  test  section  is  14  inches  in  diameter.  The  lengths 
of  ducting  were  designed  in  accordance  with  the  ASME  Power  Test  Code 
[7]  for  establishment  of  uniform  flow  and  are  shown  in  Fig.  4. 

The  entrance  to  the  test  duct  is  on  the  exterior  of  the  building 
to  minimize  thermal  gradients  in  the  air  stream.  The  entrance  is  three 
feet  square  and  is  covered  by  a  fine  mesh  screen  to  eliminate  foreign 
matter  from  the  air  stream.  The  transition  from  three  feet  square  to 
12  inches  square  is  made  in  the  entrance  section,  each  side  of  which  has 
the  curvature  of  a  quarter  ellipse.  After  a  five  foot  length  of  duct, 
a  second  transition  Is  made  from  12  inches  square  to  the  test  core 
dimensions.  This  is  followed  by  the  forward  instrument  section  which 
contains  four  thermocouple  taps,  a  piezometer  ring  and  two  pitot  tube 
taps  for  conducting  vertical  and  horizontal  velocity  surveys.  Follow¬ 
ing  the  test  core  is  the  after  instrument  section  which  contains  nine 
thermocouple  taps  and  a  piezometer  ring.  A  transition  is  then  made 
a  14  inch  outside  diameter  duct  at  the  end  of  which  is  a  standard  ASME 
square  edged  orifice.  Piezometer  rings  are  located  one  diameter 
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upstream  of  the  orifice  and  one-half  diameter  downstream  of  the  orifice. 
The  duct  section  downstream  of  the  orifice  also  contains  a  thermocouple 
tap. 

Air  flow  is  regulated  by  a  manually  operated  double  sliding  plate 
valve  located  at  the  turbo-compressor  inlet.  The  valve  plates  slide  on 
Teflon  runners  and  the  drive  screw  is  ball  bearing  mounted  to  eliminate 
maintenance  problems.  The  double  plates  constrict  the  flow  symmetrically 
and  facilitate  fine  control  of  the  air  flow.  Coarse  control  of  the  air 
flow  is  accomplished  with  the  blast  gate  (discharge  valve)  of  the  cen¬ 
trifugal  compressor. 

Air  is  induced  through  the  test  section  ducting  by  a  centrifugal 
compressor  rated  6000  cfm  at  2.5  psi  back  pressure.  The  compressor  is 
driven  by  a  three  pvsse,  220  volt,  100  hp  motor.  The  compressor  dis¬ 
charges  air  through  a  20  inch  square  duct  to  the  exterior  of  the 
building. 

The  air  flow  is  metered  by  means  of  a  standard  ASME  square  edged 
orifice  designed  to  the  specifications  given  in  the  ASME  Power  Test 
Code  t'7j .  Four  orifice  plates  made  of  1/4  inch  type  304  stainless  steel 
were  designed  with  the  following  orifice  diameters.  The  preceeding  duct 
has  an  inside  diameter,  D,  of  13.875  inches. 


d  (in) 

a  m  d 
8  -  Tf 

10.406 

0.75 

6.244 

0.45 

3.468 

0.25 

2.081 

0.15 

These  diameters  were  selected  so  that  flow  rates  of  from  250  lbs  per  hour 
to  22,000  lbs  per  hour  can  be  measured  with  the  ranges  overlapping.  The 
orifice  plates  are  fitted  with  a  centering  grove  on  the  bottom  and  a 
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scribe  mark  on  the  top  so  that  they  may  be  accurately  positioned. 

The  ducting  downstream  of  the  test  section  is  insulated  with  two 
inch  fiber  glass  insulation  covered  by  aluminum  foil,  and  the  air  duct 
is  insulkted  from  the  test  core  by  1/8  inch  Teflon  gaskets. 

Steam  System.  Heating  of  the  test  core  is  accomplished  by  passing 
saturated  steam  through  one  side  of  the  core.  Before  entering  the  core 
the  steam  is  reduced  to  a  low  pressure,  about  six  psig,  and  a  slightly 
superheated  state,  about  five  degrees  superheat.  In  a  superheated 
condition,  the  state  of  the  steam  can  accurately  be  established  by 
pressure  and  temperature.  The  few  degrees  of  superheat  are  quickly 
lost  in  the  test  core  and  have  negligible  effect  on  the  calculations. 

A  considerable  excess  of  steam  is  passed  through  the  core  to  prevent 
a  thick  film  of  cpndensate  from  building  up  on  the  heat  transfer  sur¬ 
faces. 

Steam  is  supplied  from  the  central  steam  heating  system. at  125 
psig  in  a  saturated  state.  It  is  strained  and  then  reduced  to  about 
30  psig  through  a  1  -  1/A  inch  air  operated,  pilot  controlled  pressure 
reducer.  At  this  point  water  is  injected,  if  necessary,  to  desuperheat 
the  steam  and  then  it  flows  through  a  1/4  inch  centrifugal  separator 
which  removes  all  entrained  moisture.  The  steam  then  passes  through  a 
second  ;1  -  1/4  inch  air  operated,  pilot  controlled  pressure  reducer 
which  reduces  it  to  about  six  psig.  This  second  reduction  slightly 
superheats  the  steam.  Just  prior  to  entry  into  the  core  the  steam 
passes  through  a  short  transition  section.  To  provide  uniformity  of 
flow  in  the  test  core,  the  transition  section  is  filled  with  aluminum 
shavings  held  in  place  by  a  fine  mesh  stainless  steel  screen.  Steam  and 
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condensate  from  the  bottom  of  the  test  core  flow  through  a  second 
transition  section  to  a  1  -s  1/4  Inch  centrifugal  separator  where  the 
condensate  is  removed.  A  float  type  trap  receives  the  condensate  and 
discharges  it  at  a  realitively  constant  flow  rate.  The  cbndensate  is 
then  subcooled  in  a  small  tap  water  counter  flow  heat  exchanger  and 
collected  in  a  bucket  for  weighing.  The  exhaust  steam  flows  from  the 
separator  through  a  standard  ASME  square  edged  orifice  with  flange 
pressure  taps  and  is  discharged  to  the  atmosphere. 

Two  orifice  plates  were  selected  to  measure  the  steam  flow  rates 
and  have  the  following  orifice  diameters.  The  inside  pipe  diameter, 
Ds,  is  equal  to  1.25  inches. 


0.700  .560 

0.890  .712 

The  metering  orifice  was  installed  as  specified  in  the  ASME  Power 
Test  Code  {7j. 

Pressure  and  thermocouple  tap3  were  installed  in  both  the  upper 
and  lower  transition  sections  and  a  thermocouple  tap  was  installed 
upstream  of  the  steam  orifice.  The  entire  steam  system  is  well  insur- 
lated  to  minimize  heat  loss. 

Compressed  air  for  actuation  of  the  two  steam  pressure  reducers 
is  supplied  by  two  air  pressure  regulators  located  on  the  instrument 
panel.  This  combination  is  designed  to  hold  the  downstream  steam 
pressure  to  it  0.2  psig.  An  air  line  ftom  the  air  supply  of  the  first 
pressure  reducer  was  fitted  into  the  steam  strainer  clean-out  plug  so 
that  the  entire  steam  system  can  be  blown  dry  upon  completion  of  a 


13 


series  of  test  runs.  This  should  minimize  rusting  when  the  system  stands 
Idle  for  long  periods  of  time. 
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Pressure  Instrumentation.  Pressure  instrumentation  is  provided 
for  measuring  ,  the  gage  pressure  in  the  air  duct  upstream  of  the  test 
core  and  upstream  of  the  air  orifice,  the  air  pressure  differential 
across  the  test  core  and  air  orifice,  the  steam  header  pressures  and 
the  pressure  differential  across  the  steam  orifice.  Air  pressures  are 
measured  with  well-type  single  leg  water  manometers  and  inclined  draft 
gages.  Steam  pressures  are  measured  with  well  type  single  leg  mercury 
manometers. 

A  one  inch  and  three  inch  inclined  draft  gage  and  a  60  inch  water 
manometer  are  ised  for  measuring  the  air  pressure  differential  across 
the  test  core.  A  30  inch  water  manometer  is  used  for  measuring  gage 
pressure  upstream  of  the  test  core.  The  pressure  tap  downstream  of  the 
test  core  is  located  so  that  lull  pressure  recovery  of  the  air  stream 
is  attained.  The  air  orifice  pressure  differential  is  measured  with  a 
three  inch  inclined  draft  gage  and  a  30  inch  water  manometer.  A  60  inch 
water  manometer  is  used  to  measure  the  gage  pressure  upstream  of  the 
air  orifice.  Each  pressure  tap  os  the  air  duct  consists  of  four  1/16 
inch  holes  drilled  symmetrically  around  the  duct  connected  together  by 
a  piezometer  ring  of  1/4  inch  copper  tubing.  One  quarter  inch  plastic 
tubing  is  used  to  connect  the  piezometer  rings  to  -he  manometers  and 
draft  gages  and  each  manometer  and  draft  gage  is  isolated  by  a  brass 
toggle  valve. 

Three  30  inch  mercury  manometers  are  used  to  measure  the  gage 
pressure  of  the  top  and  bottom  steam  headers  and  the  pressure  differential 
across  the  steam  orifice.  Water  is  used  as  the  intermediate  fluid 
between  the  steam  and  mercury  and  a  constant  head  of  water  is  maintained 
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by  using  water  pots.  The  change  of  level  in  the  water  pots  is  the  same 
order  of  magnitude  as  the  change  in  level  in  the  mercury  wells  and  these 
cnanges  have  insignificant  effect  on  the  pressure  readings.  The  water 
pots  are  stainless  steel  cylinders  three  inches  in  diameter  by  three 
inches  high  and  are  located  above  the  top  of  the  manometers  so  that  air 

i  < 

pockets  will  not ‘form  in  the  water  lines  connecting  them  to  the  manometers, 

\ 

The  connecting' water  lines  are  of  1/4  inch  stainless  steel  tvoing  and 
each  manometer  is  isolated  by  a  stainless  steel  toggle  valve.  The .water 
pots  are  connected  to  the  steam  headers  and  steam  orifice  by  5/8  inch 
copper  tubing.  The  large  diameter  tubing  was  selected  so  that  ecu  den¬ 
tate  can  flow  freely  back  into  the  steam  system  and  not  form  pockets 
which  would  change  the  static  head  on  the  manometers. 

Temperature  Instrumentation.  All  temperatures  are  measured  with 
copper-cons tantan  thermocouples.  The  thermocouples  are  in  assemblies 
manufactured  by  Honeywell  under  the  trade  name  MegapaK.  The  MegapaK 
consists  of  a  "head"  for  connecting  the  extension  wires,  a  "sheath"  of 
1/8  inch  stainless  steel  tubing  with  the  insulated  leads  extending  down 
through  the  center  and  the  measuring  junction  at  the  end  of  the  sheath. 

The  thermocouples  used  for  measurir„  the  air  stream  temperatures  ar« 
"exposed,"  meaning  that  the  junction  is  extended  one  sheath  diameter 
past  the  end  of  the  sheath.  Those  used  for  measuring  the  steam  tempera¬ 
tures  are  "remote,"  meaning  that  the  junction  is  one  sheath  diameter 
from  the  end  of  the  sheath  and  the  end  of  the  sheat  is  sealed  for  pressure 
tightness.  All  extension  wire  is  polyvinyl  covered  24  gage  copper-con- 
s tantan. 

The  thermocouple  voltages  are  measured  with  a  Honeywell  24  point 
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Multipoint  recorder.  The  reference  junction  is  an  ice  junction  and  is 
included  in  the  measuring  circuit.  The  recorder  is  a  D.C.  millivolt 
recorder  with  a  scale  of  zero  to  ten  millivolts.  The  smallest  division 
on  the  chart  paper  is  1/10  of  a  millivolt  and  temperatures  may  be 
estimated  to  1/100  of  a  millivolt,  The  recorder  has  several  modes  of 
operation.  On  the  Multipoint  setting  it  will  balance  and  print  each  of 
the  24  points  consecutively  at  a  rate  of  about  two  seconds  per  point. 

On  the  SelectorPrint  setting  the  recorder  will  balance  and  print  at 
about  two  second  intervals  any  number  of  the  24  points.  On  the  Hold-on 
seating  it  continuously  balances  on  any  one  of  the  24  points  but  will 
not  print.  This  setting  may  be  used  for  obtaining  a  pen  trace  by 
attaching  a  clip-on  pen  to  the  print  mechanism.  The  recorder  should 
be  adjusted  for  zero  and  span  before  use  ana  the  adlustment  procedure 
is  given  in  the  operating  procedures. 

The  temperatures  of  the  air  upstream  and  downstream  of  the  core, 
and  the  downstream  of  the  orifice  are  recorded.  Steam  temperature  in 
the  top  and  bottom  headers  and  upstream  of  the  steam  orifice  are  also 
recorded.  To  obtain  an  average  temperature  of  air  upstream  and  down¬ 
stream  of  the  core  a  system  of  thermocouples  is  used.  Upstream  of  the 
core  this  consists  of  three  thermocouples  placed  on  the  horizontal  center 
line.  Downstream  of  the  core  nine  thermocouples  are  arranged  in  a  three 
by  three  grid  with  e4ual  spacing  between  each  of  the  thermocouples  and 
the  thermocouples  and  the  wall.  In  this  manner  small  variations  in 
temperature  across  the  stream  can  be  averaged  out.  All  other  tempera¬ 
tures  are  measured  by  a  single  thermocouple.  Each  penetration  through 
the  duct  for  a  thermocouple  is  fitted  with  a  Swagelok  compression  fitting 


so  Chat  Che  penetration  length  can  be  adjusted. 

The  thermocouples  are  numbered. as  shown  in  Figs.  4  and  5  and  the 
same  number  is  printed  beside  the  temperature  reading  on  the  multipoint 
recorder  chart. 


4.  Experimental  Method  and  Equations 

General .  The  experimental  method  and  equations  used  are  essentially 
those  outlined  by  Kays  [3]  for  a  similar  facility.  These  equations  are 
"*t>4££d  for  use  with  tables  and  graphs.  The  equations  outlined  herein  are 
suited  for  use  in  a  digital  computer  program  which  will  accept  raw  test 
data  and  produce  finished  results.  The  equations  used  la  preparation  of 
original  graphs  and  tables  were  used  if  available  and  approximate 
equations  were  written  for  those  not  availably.  Graphs  of  the  approxi¬ 
mations  are  shown  in  the  appendix  and  test  data  should  be  checked  for 
applicability  of  the  approximations. 

Mass  Rate  of  Air  Flow.  The  mass  rate  of  air  flow,  m,  through  an 
ASME  square  edged  orifice,  is  given  in  the  Power  Test  Code  {7j  as: 

m  =  359  C  F  d2  F^Y  V$P0f0  '  lbm/hr 


where 


C 

F 

d 

FA 

Y 

APC 

Po 


*  coefficient  of  discharge 

■  velocity  of  approach  factor 

■  orifice  diameter,  in. 

■  orifice  thermal  expansion  factor 

*  net  expansion  factof 

■  pressure  differential  across  the  orifice,  in.  I^O 

«  specific  weight  of  the  air  upstream  of  the  orifice, 
lbp/ft3 
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The  ratio  of  the  orifice  diameter  to  the  duct  diameter  is: 

0  -  i 

D 

where  D  »  duct  diameter 

-  13.875  in. 

and  the  velocity  of  approach  factor  is: 

F  ■ 

The  absolute  duct  pressure  upstream  of  the  orifice  xs  given  by: 

po  =  pbC0,490)  -  P0  (0.0361)  psid 

where  P^  -  barometric  pressure,  in.  Hg 

Pc  -  duct  pressure  upstream  of  the  orifice,  in.  1^0 
(duct  pressure  is  below  atmospheric  pressure) 

A  dimensionless  ratio  of  orifice  pressure  drop  to  orifice  pressure 

is: 

x  » 

and  the  net  expansion  factor  is  then: 

Y  -  1.0  -  4-  0.3  5  p4) 

The  approximate  equation  for  the  orifice  thermal  expansion  factor 
for  type  304  stainless  steel  gives: 

FA-  1.0  4- (to-  68)0-85  xio"5) 

The  coefficient  of  discharge  is  an  unknown  in  the  equation  but  may 
be  quickly  iterated  using  an  equation  suggested  by  Murdock  [6].  The 
equation  is  dependant  upon  the  Reynolds  number  in  the  orifice,  NRe(j, 

/  ».r\4 


and 
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The  coefficients  CQ  and  AC  for  applicable  6's  from  Murdock  [6]  and  a 
suggested  first  iteration  C  are: 


0 

CD  AC 

C 

.15 

0.59446 

0.00945 

0.59/75 

.25 

0.59483 

0.01037 

0.5966 

.45 

0.59863 

0.01543 

0.6014 

.75 

0.60480 

0.05448 

0.6128 

The  density  of  the  air  upstream  of  the  orifice  may  tye  found  from 
the  equation  of  state  of  a  perfect  gas  modified  by  a  humidity  correction 


factor  which  gives: 


fc- 

where  P0  * 
To  - 


144  Xpn 
53.3  4  T0 

duct  pressure  upstream  of  the  orifice,  psia 

temperature  of  air  upstream  of  the  orifice,  °R 

humidity  correction  factor  for  density  as  given  by 
Kays  and  London  [4] 


0±hL_ 

(l  +  1.607  H) 


And  H  -  humidity  ratio,  lfo.^  Hz0/lbm  dry  air 
The  dynamic  viscosity, /Aot  for  air  at  one  atmosphere  pressure  is 
given  in  reference  (1)  as: 


}Xo  - 


3/2 

-  0*00  3527  Tok 


Tok  +  110.4 


where  T0^  ■  temperature  of  air  in  duct,  °K 
The  orifice  Reynolds  number:  may  be  found  from: 


N 


Red  ~ 


15,2  8  ft 

Mo  * 


Heat  Transfer  Calculations.  The  unit  conductance  for  thermal  con¬ 


vection  heat  transfer  as  given  by  equation  one  in  the  theory  section  is: 

1 
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where 
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Aa 
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Ac 


*v 


unit  conductance  for  thermal  convection  heat 
transfer  to  the  air,  BTU/ (hr  ft2  °F) 

total  surface  effectiveness  for  the  air  side 

overall  unit  conductance,  BTU/ (hr  ft2  °F) 

total  heat  transfer  area  on  the  air  side,  ft2 

total  surface  effectiveness  for  the  steam  side 

total  heat  transfer  area,  on  the  .steam  side,  ft* 

unit  conduction  for  thermal  convection  on  the 
steam  side,  BTU/(hr  ft2  op) 

thickness  of  the  wall  separting  the  steam  side 
from  the  air  side,  in. 


Awa  »  area  of  the  wall  on  the  air  side,  ft/ 

*  thermal  conductivity  of  the  wall,  BTU/ (hr  ft2  °F/ft) 
The  film  coefficient  for  convection  on  the  steam  side,  hg,  is  unknown 
and  must  be  estimated.  Any  reasonable  estimate  is  considered  adequate 
since  a  50  percent  error  in  the  estimate  results  in  about  two  percent 
error  in  ha.  A  value  of  hs  equal  to  2000  BTU/ (hr  ft2  °F)  has  been  used. 

The  specific  heat  of  dry  air  may  be  considered  constant  and  is  assumed 
to  be: 


Cpcte  =  0.241  BTU/Obr,,  *F) 


The  humidity  correction  factor  for  the  specific  heat  as  given  by  Kays 
and  London  [4]  is: 

1  +  1*91 5H 

Xc“  I  +  H 

where  H  ■  humidity  ratio,  lb  m  ^O/lb^  dry  air 
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and  the  corrected  specific  heat  is  then: 

Cp  =  Xc  BTl^tik^  F) 

The  bulk  average  temperature  in  the  test  core  is: 


whp 


re 


t>b  = 
ti  - 
t2  - 


"fci  1"  #r 

i— . — -u—  F 


f\ 

c. 


air  temperature  upstream  of  the  test  core,  °F 
air  temperature  downstream  of  the  test  core,  °F 


and  the  film  average  temperature  is  then: 

o>- 


tf=  ^>4^ 


where  ts  -  the  saturation  temperature  corresponding  to  the 

arithmetic  average  steam  pressure  in  the  core,  °F 

The  number  of  heat  Transfer  Units  may  be  found  from  the  temperature 

NTU  relation  as  given  in  the  theory  section,  yielding: 

r  te¬ 


rn 


^s*~t|  \ 

t2/ 


and  the  overall  .unit  conductance  is  then: 

y  -  cPmNTU  BTu/(hr  ft2  CF) 

For  most  plate-fin  type  geometries  the  equation  for  the  temperature 
effectiveness  of  a  straight  fin  with  constant  conduction  cross  section 
may  be  used  to  a  good  approximation. [3] .  The  fin  temperature  effective¬ 
ness  for  the  steam  side,  rif 8 ,  is: 


'W  - 


tanh(nn$*s,/l,2) 


where 


effective  fin  length  fox  the  steam  side  and  Is 
equal  to  one-half  the  total  fin  height,  in. 


V 


24  hs 
kfs 


ft 


-I 
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and  hc  *  2000  BTU/(hr  ft^  °F)  as  previously  estimated. 

kgs  =  thermal  conductance  for  the  fin  on  the  steam 
side,  BTU/ (hr  ft2  °F/f t) . 

6g  ■  fin  thickness  on  the  steam  side,  in. 

The  overall  temperature  effectiveness  for  the  steam  side,  ns,  is  then: 


o 

where  A^s  *  wall  area  on  the  steam  side,  ftA 

2 

Afg  *  fin  area  on  the  steam  side,  ft 

2 

As  *  total  heat  transfer  area  on  the  steam  side,  ft 
The  fin  temperature  effectiveness  for  the  air  side,7||‘a,  is; 

tanh  (ma*a/12) 

where  aa  •=■  effective  fin  length  for  the  air  side  and  is  equal 

to  one-half  the  total  fin  height,  in. 


and  kfs  =  thermal  conductance  for  the  fin  on  the  air  side, 

BTU/ (hr  ft2  °F/ft). 

oa  «*  fin  thickness  on  the  air  side,  in. 

ha  »  film  coefficient  for  convection  on  the  air  side. 

It  may  be  taken  equal  to  U  for  the  first  iteration 
BTU/ (hr  ft2  oF) , 

The  overall  temperature  effectiveness  for  the  air  side,  nQ>  is  given  by: 


where  Awa  -  wall  area  on  the  air  side,  ft2 
Af4  *  fin  area  on  the  air  side,  ft2 
Aa  »  total  heat  transfer  area  on  the  air  side,  ft2 
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The  air  mass  velocity  in  the  core,  G,  is  given  by: 


o  .  A  lbm//hr  ft2) 

Ac 

where  Ac  *  the  free  flow  area  of  the  core  on  the  air  side,  ft^ 
Stanton  number  in  the  core  can  now  be  evaluated  as: 

NSt  =  7^~ 

5t  Scp 

The  dynamic  viscosity  evaluated  at  the  film  average  temperature  and 
one  atmosphere  pressure  as  given  by  Hilsenrath  [1]  is: 

3/2 

u  .  0.00?$27  %  it  /(hr  ft) 

Tfk+-  110.4  m/' 


where  T^  -  film  average  temperature,  °K 

The  dynamic  viscosity  evaluated  at  the  bulk  average  temperature  is: 

3/> 

0.003527  Tuf  ft) 

/^b  Tkl.  +■  110.4  w  '  1 


where  Tbk  ■  bulk  average  temperatupe,  °K 

Ther  thermal  conductivity  of  air  at  the  film  average  temperature 
and  one  atmosphere  pressure  as  set  forth  by  Hilsenrath  [1]  is: 

T 

.  ,  0.001529  Tfk  BTU _ 

f  r,+ .  2<s-<  .  i  " ft*  F/ft 


245.4 


Where  T^  ■  film  average  temperature,  °K 

The  Reynolds  numbers  evaluated  at  film  average  temperature  and  bulk 
average  temperature  are  evaluated  by: 

_  4rhG 


Ref  "  JUf 
'Reb  "  yub 
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where  hydraulic  radius  rjj  is  defined  by  Kays  and  London  [3]  to  be  based 


on  the  free  flow  area  and  takes  the  form 


where  L  «  core  flow  length,  ft 
Prantdl  number  in  the  core  is: 


N 


_  Mfc 

Pr  kf 


E 


where  the  fluid  properties  are  evaluated  at  the  film  average  temperature 
in  the  test  core. 

The  Colburn-j  heat  transfer  modulus  is  then: 


N 


Pr 


which  concludes  the  heat  transfer  calculations. 


Friction  Factor  Calculations.  The  equation  for  the  core  friction 
factor  as  developed  in  equation  two  in  the  theory  section  may  be  written 


as’  f  rh  Pm  . 

y  APC(4.3255  X  I09) 

Kc  -("I  +  (T2) 

L 

f  G2 

u 

L  ?i  J 

where  r^ 
Pm 
L 

APC 

6 


core  hydraulic  radius,  ft 

mean  density  in  the  core  lbm/ft^ 

core  flow  length,  ft 

pressure  drop  across  core,  in.  1^0 

core  mass  velocity,  lbm/hr  ft^ 
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contraction  coefficient 


Kc  * 

CT~  “  the  ratio  of  free  flow  area  to  frontal-area 
*  density  of  the  air  upstream  of  core,  lb^ft^ 

K%  ■  expansion  coefficient 

^2  *  density  of  the  air  downstream  of  core,  lbm/ft^ 

fd  *  friction  factor  for  the  duct  downstream  of  the 
core 

1^  *  length  of  duct  from  the  core  to  the  downstream 

pressure  tap,  in. 

D(j  *  hydraulic  diameter  of  the  duct  downstream  of 
core,  in. 

The  duct  friction  factor,  fd,  is  assumed  constant  at  .015,  a  value 
correspondingy.to  a  nearly  smooth  pipe  for  a  Reyholds  number  range  of 
100,000  to  500,000  as  may  be  expected  in  the  downstream  duct. 


The  mean  density  in  the  case  of  heating  is  given  by: 

T,e, +T2f2 


2 


upstream  temperature,  °R 
downstream  temperature,  °R 

Saturation  temperature  corresponding  to  the  average 
steam  pressure  in  the  core  °R 


upstream  density 


° 2  ■  downstream  pressure,  in.  H2O  absolute 

In  the  isothermal  case,  the  mean  density  is  equal  to  the  arithmetic 
average  density  and  since  the  temperature  remains  constant 

.  0.097318  lb  /ft8 

-  Cave-  2Tj 

The  equations  for  the  contraction  coefficient,  Kc,  and  the  expansion 
coefficient.  Kg,  are  developed  by  Kays  [2]  and  are  given  by 


Kr -- 


I  -  2C--+-  C, 


(gXd-0 


Ke  =  1  -  <r  +-  <r‘ 


where  *  velocity-distribution  coefficient 

<r  -  Ac/Afr 

Cc  =  jet  contraction  ratio 

The  jet  contraction  ratio  is  approximated  as 

Cc  =  0.611  4-0.45  <r+  0.344  (T5,7 

The  velocity-distribution  coefficient  K<j  is  dependent  upon  both  core 


geometry  and  film  Reynolds,  NRef  The  usual  geometries  used  in  plate 


fin  type  construction  are: 


square 

triangular 


For  laminar  flow,  %ef  less  than  2000,  Kays  1 2 ]  gives  the  values 


square  Kd  ■  Kj  -1.39 

U  Ug 

triangular  Kjj  -  K<jt  -  1.43 

and  for  turbulent  flow 
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square  Kj 


1  +  1.17  (Kdc  -  1) 


triangular  K^  «  K<jt  *  1  +  1.29  (K^  -  1) 

where  the  velocity-distribution  coefficient  for  circular  tubes,  X,jg,  is 


given  by: 


Kdc-  U09068  (4fm)  +  0.0584  ^4fJ~fT 


fm  *  0.049  NRef 


-0.2 


If  the  boundary  layer  is  continually  tripped,  as  in  the  case  of  louvered 
fin  or  off-set  fin  construction,  film  Reynolds  number  should  be  taken 
equal  to  infinity  which  yields  Kd  *»  1.0  in  all  cases. 


Lculations.  The  heat  transfer  rate  to  the  air. 


Q  air,  is: 


Q  air  ■=  m  Cp  (t2  -  ti)  BTU/hr 


The  total  heat  given  up  by  the  steam,  Q  steam,  is  equal  to  the  heat 

given  up  by  the  excess  steam  as  it  passes  through  the  core  plus  the  heat 

given  up  by  the  condensate  in  the  core. 

Q  steam  =  s  (hi  **  h2g)  +  wc  (hi  -  l^f)  BTU/hr 

where  s  **  mass  rate  of  flow  of  excess  steam,  lbm/hr 

wc  »  mass  rate  of  flow  of  the  condensate,  lbjjj/hr 

h->  *  enthalpy  of  the  steam  at  the  top  header  pressure  and 
temperature  BTU/lbm 

enthalpy  of  saturated  steam  at  the  bottom  header 
pressure,  BTU/lbm 

h2f*  enthalpy  of  saturated  liquid  at  the  bottom  header 
pressure,  BTU/lbm 

The  mass  rate  of  flow  of  the  excess  steam,  s,  as  measured  by  a  square 
edged,  standard  ASMS  orifice  with  flange  pressure  taps  is  found  from  the 
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applicable  equation  given  in  the  ASME  Power  Test  Code  [7}: 


S  =  359  C*  Fs  df  FAsYs  ^  Ps  '  lbm/  hr 

where  Cg  ■  coefficient  of  discharge 

Fs  *  velocity  of  approach  factor 
de  «*  orifice  diameter,  in». 

Fas  *  orifice  thermal  expansion  factor 
Yg  *  net  e. 'pension  factor 

Aps  *  pressure  differential  across  the  s  .  »am  orifice, 
in.  H20 

Vg  *»  specific  volume  of  the  steam  passing  through  the 
orifice,  ft^/lbm 

Tne  velocity  of  approach  factor,  Fs,  is: 


where 


and  Ds  *  steam  pipe  diameter 

«  1.25  in. 

The  orifice  thermal  expansion  factor  is  taken  as  constant  corres¬ 
ponding  to  a  steam  temperature  of  220  °F.  For  type  304  stainless  steel 


1.00.'/ 


The  net  expansion  factor  for  steam,  Ys,  is: 

YS~|-(0.4I  -FO«35054)rxs/l*3) 

where  '  ' 


XS®  APs/’Ps 

The  upstream  steam  pressure,  Pg,  is  taken  as  the  saturation  pressure 
corresponding  to  the  upstream  temperature  and  is  given  by  the  approximate 
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equation 


P*=  16.5  +  (tso  “218X0.376)  psta 

The  specific  volume  is  given  by  the  approximate  equation: 

Vg  =  24 #01  0.4 (tso  2 18)  ft/ lb 

The  coefficient  of  discharge,  C  ,  is  iterated  by  the  same  method 
that  was  used  for  the  air  orifice.  The  values  of  CQ  and  AC  for  applicable 
values  of  $s  are: 


^s 

'  ds 

Co  "  Cos 

AC  -  AC 

.56 

.700  in. 

.6085 

.025 

.71 

.890  in. 

.6090 

.041 

and  then 

Cs=-Cos't" 

(m°4  f 

The  steam  orifice  Reynolds  number,  NRe8(j,  is: 


/“S  els 


and  the  dynamic  viscosity  is  given  by  the  approximate  equation: 

8.25  X  KfVfco-^OOXl-^  x  to8) 


ft  sec 


The  enthalpy  of  the  superheated  steam  in  the  top  header  is  approxi¬ 
mated  by  two  equations,  the  first  giving  the  variation  with  tempqrdture 
at  a  constant  pressure,  P8^  ■  18.0  psia  and  the  second  equation  the 
variation  with  pressure.  These  equations  are: 

h^  -  1158.0  +  (t8l  -  230.0)  (.5)  BTU/lbm 

where  h^  «  enthalpy  at  13.0  psia> 

t8l  *  top  header  steam  temperatrue,  °F 


35 


-  *  >-'*%■**!***•'.-« ■*&**+  -i,  mewaata ¥V**amo»tfmm  wwBWatww«Wlft 


thence  hsl  -  h/j.  -  (PsjL  -  18.0)  (.35)  BTU/lbn 

where  PSX  *  steam  pressure  in  top  header,  peia 

The  enthalpy  of  the  saturated  st^am  in  the  bottom  header  is  found 
from  the  approximate  equation: 

h82  “  1155.3  +  (Po2  -  19.0)  (.934)  BTU/lbn 

where  Ps2  *  stream  pressure  in  bottom  header,  psia 

and  the  enthalpy  of  the  condensate  in  the  bottom  header  is: 

hc  -  193.4  +  (Ps2  -  19.0)  (2.64)  BTU/lbm 

The  mass  rate  of  flow  of  the  condensate,  wc,  is  found  by  timing 
the  condensate  collected. 

5.  Evaluation 

A  steady  state  steam  to  air  test  was  conducted  on  a  compact  heat 
exchanger  manufactured  by  Hamilton  Standard  Corporation.  The  test  core 
had  the  approximate  dimensions  of  a  six  inch  cube.  The  construction  was 
the  plane  plate-fin  type  with  flow  passages  between  fins  being  approxi¬ 
mately  square.  One  side  of  the  core  was  '’straight,"  meaning  that  the  flow 
passage  was  straight  for  the  length  of  the  exchanger.  The  other  side  was 
"ruffled,"  meaning  that  the  flow  passage  was  approximately  a  cosine  curve. 
The  plate  spacing  was  0.156  inches  and  the  fin  spacing  was  18  fins  per 
inch  for  both  sides.  All  plate  and  fin  material  was  0.004  inch  stainless 
steel.  A  set  of  "standard"  curves  for  Colburn-j  and  friction  factor  vs. 
Reynolds  number  were  supplied  with  the  test  core.  The  "standard"  curves 
were  for  plate  and  fin  spacing  similar  to  the  core  tested. 

The  test  was  performed  using  the  straight  side  of  the  core  for  air 
and  the  ruffled  side  for  steam.  The  various  components  of  the  facility, 
the  air  system,  steam  system  and  instrumentation  were  checked  out  and 
proved  to  be  very  satisfactory.  There  was  a  problem  in  defining  flow  and 
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heat-  transfer  areas  for  the  exchanger  because  structural  stiffeners  had 
been  placed  at  the  end  of  each  plate-fin  passage.  The  stiffeners  were 
placed  at  an  angle  across  the  flow  path  and  thereby  restricted  the 
flow  in  the  peripheral  flow  passages.  This  restriction  w..uld  tend  to 
increase  with  increasing  flow  rates.  The  test  results  based  upon  the 
minimum  free  flow  area,  which  neglects  the  peripheral  flow  passages 
affected  by  the  stiffeners,  was  plotted  on  the  standard  curves  as 
shown  in  Fig.  6. 

The  next  test  should  be  on  a  heat  exchanger  core  with  previously 
determined  characteristics  and  with  given  or  definable  geometry. 
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APPENDIX  I 
OPERATING  PROCEDURE 

General ,  The  preliminary  check  of  instrumentation  for  both  the 
hot  core  and  cold  core  tests  is  the  same. 

The  multipoint  recorder  should  be  allowed  to  warm  up  two  hours 
prior  to  calibration.  To  calibrate  the  instrument,  set  the  print 
mechanism  to  "Hold  On"  indicating  any  one  of  the  points  not  used  for  a 
thermocouple.  This  open  circuits  the  thermocouple  measuring  circuit. 
Place  a  D.C.  millivolt  source  across  the  terminals  at  the  back  of  the 
recorder  and  set  the  source  to  one  millivolt.  Adjust  the  recorder  zero 
so  that  one  millivolt  is  indicated.  Set  the  source  to  nine  millivolts 
and  adjust  the  span  so  that  nine  millivolts  are  indicated.  Repeat  the 
zero  and  span  adjustments  until  the  recorder  is  calibrated. 

All  manometers  should  be  adjusted  for  zero  with  the  isolation  valves 
open  prior  to  lighting  off  the  air  or  steam  systems.  The  two  steam 
header  ir->rcury  manometers  must  be  zeroed  at  about  three  Inches  due  to  the 
head  of  water  between  the  manometers  and  the  water  pots.  The  water  lines 
between  the  mercury  manometers  and  the  water  pots  are  susceptible  to  the 
formation  of  air  pockets  which  will  cause  erroneous  readings. 


Cold  Core  Test.  In  order  to  assure  an  even  distribution  of  data 


points,  the  approximate  pressure  differential  across  the  air  orifice 
should  be  calculated  for  each  run.  Calculate  the  mass  rate  of  air 
flow  at  a  core  Reynolds  number  of  500  and  10,000  and  plot  these  two 
points  on  a  log  log  plot.  Divide  the  straight  line  connecting  the  points 
into  equal  increments,  the  number  Spending  upon  the  desired  number  of 
runs.  The  orifice  pressure  differential  may  then  be  read  from  Fig.  7 
for  each  mass  rate  of  air  flow. 

Open  all  manometer  valves  and  adjust  for  zero.  Close  all  manometer 
valves  prior  to  lighting  off.  Set  the  multipoint  recorder  to  print  the 


thermocouples  in  use  by  raising  the  correapt/nuxii)^  uullous  on  oexecto** 


Print."  Turn  the  chart  motor  on  and  select  the  desired  speed.  Close 
the  compressor  blast  gate  (discharge  valve)  and  energize  the  compressor 
motor.  Open  the  blast  gate  to  desired  setting.  Adjust  the  sliding 
plate  valve  for  the  desired  air  orifice  pressure  differential  and  observe 
instrumentation  until  steady  state  obtained. 

Record  the  following  data  for  each  run: 

1.  Atmospheric  pressure  (0.001  in, Hg) 

2.  Atmospheric  dry  bulb  temperature  (0.1  °F) 

3.  Atmospheric  wet  bulh  temperature  (0.1  °F) 

A,  Air  orifice  size 

5.  Orifice  upstream  pressure  (0.01  in,  1^0) 

6.  Orifice  pressure  differential  (0.001  in,  H20  for  APo  <3",  and 
0.01  in.  H20  for  APg  >3") 

7.  Core  upstream  pressure  (0.01  in,  H20) 

8.  Core  pressure  differential  (0.001  in, H20  for  APc  <3",  and  0.01 
in,  Ho0  for  APc  >3")  . 

4  O 
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9.  Orifice  air  temperature  (0.01  millivolt) 

10.  Core  upstream  air  temperature  (0.01  millivolt) 

Items  five  through  ten  are  recorded  three  times  in  the  same  sequence 
and  averaged  on  the  data  sheet. 

Hot  Core  Test.  The  instrumentation  checks  and  the  procedures  for 
lighting  off  the  air  system  are  the  same  as  for  the  Cold  Core  Test. 

The  steam  system  is  activated  by  the  following  procedures.  Obtain 
100  psig  low  pressure  air  to  the  reducer  central  valves.  Turn  on  the 
tap  water  supply  for  the  condensate  cooler  and  the  desuperheater.  Crack 
the  desuperheater  supply  valve.  Open  the  by-pass  valves  for  the  steam  j 

I 

reducers  and  crack  the  steam  stop  valve.  Allow  the  system  to  warm  up 

slowly.  When  warm-up  is  completed  close  the  by-pass  valves  and  open  the  i 

steam  stop.  Charge  the  first  steam  reducer  to  about  30  psig  and  the  ! 

i 

second  to  about  six  psig.  Obtain  a  steam  state  in  the  top  steam  header 

of  about  six  psig  with  five  to  ten  degrees  superheat  by  adjusting  the  ! 

first  steam  reducer  and  the  desuperheater.  The  excess  steam  flow  rate 

should  be  at  least  five  times  the  condensate  flow  rate.  Two  steam 

orifice  plates  are  provided;  the  larger  to  be  used  for  high  ait  flow  rates. 

Condensate  should  be  collected  for  at  least  five  minutes  at  high  air  flow 

rates  and  for  at  least  ten  minutes  at  low  air  flow  rates. 

Steady  state  must  be  achieved  before  each  run.  This  will  require 
about  oue  nour  prior  to  the  first  run  and  15  minutes  between  runs.  Steady 
downstream  core  temperatures  are  the  best  indication  of  steady  state. 

The  following  data  is  recorded  for  each  run: 

1.  Atmospheric  pressure  (0.001  in.  Hg) 

2.  Atmospheric  dry  bulb  temperature  (0.1  °F) 
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,>*■**“ . 

I 

) 

! 

3.  Atmospheric  wet  bulb  temperature  (0.1  of) 

4.  Air  orifice  size 

5.  Steam  orifice  size 

6.  Top  steam  header  pressure  (0.01  in.  Hg) 

I 

1  '  7.  Bottom  steam  header  pressure  {0.01  in.  Hg) 

l 

\  8.  Pressure  differential  across  the  steam  orifice  (0.01  in.  Hg 

under  H2O) 

9.  Core  upstream  pressure  (0.01  in.  HjO) 

10.  Core  pressure  differential  (0.001  in.  H2O  for  APc<  3",  and  0.01 
in.  H2O  for  APc>  3") 

|  11.  Orifice  pressure  differential  (0.001  in.  H2O  for  APo<  3",  and 

t  0.01  in.H20  for  APc,  >3") 

12.  Orifice  upstream  pressure  (0.01  in.  ^0) 

13.  Core  upstream  air  temperature  (0.01  millivolt) 

I 

j  -  14.  Core  downstream  air  temperature  (0.01  millivolt) 

|  15.  Orifice  air  temperature  (0.01  millivolt) 

r 

!  16.  Top  steam  header  temperature  (0.01  millivolt) 

17.  Bottom  steam  header  temperature  (0.01  millivolt) 

18.  Steam  orifice  temperature  (0.01  millivolt) 

|  19.  Condensate  weight  (0.1  1U) 

I 

j  20.  Coadensate  collection  time  (1.  secj 

Record  items  6  through  18  three  times  in  the  same  sequence  and 
average  on  the  data  sheet. 

The  humidity  ratio,  H,  may  be  found  from  a  standard  psychometric 
chart  and  convert  temperatures  to  °F  (0.1). 
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Figure  7.  Mass  rate  of  air  flow  vs* 

orifice  pressure  differential 


APPENDIX  II 


DIGITAL  COMPUTER  PROGRAM  FOR  DATA  REDUCTION 
The  progiam  is  designed  to  reduce  data  from  both  the  hof  core  test 
and  the  cold  core  test.  The  program  written  in  FORTRAN  60  and  is  shown 
on  pages  52  through  57  and  the  output  for  the  hot  core  test  is  shown  on 
page  58  and  that  of  the  cold  core  test  on  page  59  .  The  variables  used 
in  the  program  are  defined  in  the  program  glossary.  Variables  not  appear¬ 
ing  in  the  glossary  are  internal  to  the  program  and  are  used  to  change 
dimension  or  to  define  special  groupings  to  simplify  programing. 

The  index  used  to  reduce  data  from  the  hot  care  test  or  the  cold 
core  test  is  given  by 2 


Program  Index 
IHORC 

Hot  core  test  1 

Cold  core  test  2 

The  program  is  indexed  to  accept  the  various  fin  geometries  normally 
used  in  plate  Jin  construction  and  the  indices,  are  as  follows: 

Fin  Geometry  Index 
ISQTRIN 

Square  fin  1 

Triangular  fin  2 

Louvered  or  off-set  fin  3 

(Nr&  =  infinity) 

'The  indicAS  used  for  the  fluid  metering  orifices  are: 

Air  Orifice  Air  orifice  Index 

BETA  d  (in.)  IB 


45 


.75 


10.406 


4 


.45 

6.244 

3 

.25 

3.468 

2 

.15 

2.081 

1 

Steam  Orifice 

Steam  Orifice  Index 

BS 

ds  (in.) 

IBS 

.71 

0.890 

2 

.56 

0.700 

1 

The  indices  are  all  entered  as  fixed  point  constants.  The  number 
of  runs,  NOR,  and  the  run  number,  NR,  are  entered  as  fixed  point  con¬ 
stants.  All  fixed  point  constants  must  be  right  justified  in  their 
data  fields. 

The  constants  for  a  given  core  such  as  areas,  lengths,  etc.  are 
entered  as  input  data  only  once  and  are  the  first  three  dita  cards 
beginning  with  AC  and  ending  with  DD.  The  fixed  point  constants  for  a 
particular  series  of  runs,  NOR,  IHORC,  ISQTRIN  are  entered  on  the  fourth 
data  card.  The  fifth  through  the  seventh  data  cards  are  for  run  one, 
eighth  through  tenth  for  run  two  and  so  on.  The  first  data  card  for  a 
run  contains  the  run  fixed  point  constants,  NR,  IB,  IBS.  The  second  and 
third  data  cards  for  a  run  contain  the  run  data  beginning  with  K  and 
ending  with  TC.  When  entering  data  for  the  cold  '-ore  test,  only  the 
applicable  data  must  be  entered.  All  other  data  fields  may  be  left  blank. 

The  dimensions  of  input  data  must  be  those  given  in  the  program 
glossary  and  are  the  same  as  the  raw  data. 
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PROGRAM  GLOSSARY 


Program 

Program 

Nomenclatu: 

Symbol 

Unit 

Symbol 

AA 

ft2 

Aa 

AC 

ft2 

Ac 

ALFA 

AFA 

ft2 

Afa 

AFR 

ft2 

Afr 

AFS 

ft2 

Afs 

AKFM 

BTU/(hr  ft2  °F/f t) 

%m 

AS 

ft2 

As 

AWA 

ft2 

^a 

AWS 

ft2 

BETA 

6 

BS 

C 

C 

CC 

Cc 

CO 

Co 

COLJ 

j 

COS 

c 

os 

CP 

BTU/ (lbm  °F) 

c 

P 

CRKD 

Kdc 

CS 

cs 

DC 

AC 

DCS 

ACs 

DD 

in. 

Dd 

47 


I 

I 


1 

DL 

in. 

*d 

r 

5 

DPC 

in.  H20 

4Pq 

1 

DPO 

in.  H20 

AP 

0 

< 

DPS 

in.  Hg  under  water 

*Ps 

j 

i  . 

i 

i 

DSO 

in. 

ds 

i 

i 

ERROR 

Error 

} 

! 

1 

FA 

PA 

i 

FAS 

i 

fAs 

1 

FB 

F 

i 

t 

] 

j 

FBS 

Fs 

i 

t 

FC  > 

f 

i 

* 

FD 

fd 

FLA 

in. 

ta 

i 

* 

<\ 

FLL 

ft 

L 

1 

^  * 

FLS 

in. 

K 

>< 

i 

FKA 

BTU/ (hr  ft2op/ft) 

kfa 

i 

M 

FKS 

BTU/ (hr  ft2oF/ft) 

kf  s 

FM 

fm 

G 

lbs/ (hr  ft2) 

G 

l 

( 

H 

lb/lb 

H 

,f 

HA 

BTU/ (hr  f t2oF) 

ha 

HG 

BTU/lbm 

h 

c 

HS 

BTU/ (hr  f t2oF) 

hQ 

s 

HS1 

BTU/ lb 

m 

hsl 

* 

j 

HS1G 

BTU/lbm 

CO 

»-» 

HS2 

BTU/lbm 

hs2 

48 

IB 


IBS 

IHORC 

ISQTRIN 

NOR 

NR 

OD 

PI 

P2 

PB 

PO 

PR23 

PS1 

PS2 

PSO 

Q  AIR 

QSTM 

REB 

REFM 

REO 

RESO 

RH 

ROl 

R02 

ROAVE 


Index  for  air  orifice  ciameter 

Index  for  steam  orifice  diameter 

Index  for  Hot  or  Cold  core  test 

Index  for  fin  geometry,  square,  triangular 
or  infinite  Reynolds  number 

Number  of  runs 

Run  number 

in.  d 


in.  H20 

P1 

in.  H20 

P2 

in.  Hg 

Pb 

in.  H20 

Po 

N  2/3 
wpr 

in.  Hg 

Psl 

in.  Hg 

Ps2 

psia 

Pso 

BTU/hr 

Q  air 

BTU/hr 

Qsteam 

%eb 

%ef 

NRed 

NResd 

ft 

rh 

lb!  ta/f t3 

Pi 

lb:  ta/f  t3 

P2 

lb  uj/ft3 

Pave 

lb:  ta/f t- 
lb.  p/ft3 
lb  p/hr 


time  of  condensate 
collection 


BTU/ (hr  f ?°F) 
lb:  y (hr  ft) 
lb  m/ (hr  ft) 


lb.  m/ (hr  ft) 
lb: m/ (hr  ft) 
lb  m/ (sec  ft) 


yU  in  core  for  cold  core 
test 


vs 

wc 

WCH 

WDOT 

WK 

X 

xc 

SKD 

XM 

XMA 

XMS 

XNTU 

XS 

y 

YKE 

YNFA 

YNFS 

YNO 

YNS 

YS 

Z*C 


ftS/lb,,, 

Vs 

»■« 

Weight  of ■ condensate 
collected 

lb  m/hr 

9 

wc 

lb'  ta/h" 

0 

m 

BTU/ (hr  f t2oF/ft) 

x  i 

1 

XC  ] 

Kd  j 

Xm  ! 

i 

ft-1 

ma  1 

ffl 

ms  j 

NTU 

r 

xs 

X 


nfa 

nfs 

% 

n8 

Xs 


PROGRAM  SSHEAT 

OREAD.  l,AC,AFR,AA,AWA,AFA,AS,AWS,AFS,FLl,  fFS,TFA,TW, 
1FLA,FLS,FKA,FKS,WK,FD,DL,DD 

1  FORMAT (8F10 • 0 ) 

READ  2,N0R, IHQRC, ISQTRIN 

2  FORMAT (3110 ) 

GO  TO  (3,5),IH0RC 

3  PRINT  4 

40FORMAT <iHi  25X,32H  HEAT  TRANSFER  AND  FRICTION  DATA,/// 

1106H  RUN  MOOT  G  HUMID  T1  PI  T2  DPC  TS 

2H  8TU/  COL J  F  FACT  NRF  NR  B  ERROR;// 

3107H  LB/HR  L8/HR  FT2  lB/LB  F  PSIA  F  PSI  F  H 
4R  FT2  F  PERCENT,//) 

GO  TO  7 
5  PRINT  6 

60FORMAT(1H1  25X,25H  ISOTHERMAL  FRICTION  DATA,/// 

160H  RUN  MOOT  G  HUMID  T1  Pi  DPC  F  FACT  NR,// 

246H  LB/HR  L8/HR  FT2  LB/LB  F  PSIA  PSJ,//) 

7  READ  6, NR, IB, IBS 

8  FORMAT  (3110  > 

OREAD  9,H,0D,DS0,8ETA,BS,T0,Tl,T2#TSi,TS2,TS0, 
1PB,PS1,PS2,DPS,P1,DPC,PO,DPO,WC,TC 

9  FORMAT  (8F10.0) 

PS1A-«490*(PB+PS1) 

PS2As.490*(PB+PS2) 

PSAVEs(PSlA+PS2A)/2, 

TSs220,+(PSAVE-17»L)*2,66 

P2SP1+DPC 

TORsTO+459.7 

T0K=T0R*5./9. 

UO=,003527*TOK**i.5/(TOK  +  U0,4) 

GO  TO  (20,21,22,23  >,IB 

20  COs. 59446 
DCs. 00945 
Cls. 5975 
ODslO.406 
BETAs, 75 
GO  TO  25 

21  CO=. 59483 
DCs, 01037 
Cl=. 5966 
0D=6 . 244 
BETAs. 45 
GO  TO  25 

22  COs. 59868 
DCs, 01543 
Cl= ,6014 
0Ds3,468 
BETAs. 25 
GO  TO  25 

23  COs, .60480 
DC=. 05448 
Cls, 6128 
00=2,081 
BETAs, 15 

25  C=Cl 

B4sBETA**4 

POAsPB*,490-PO*,0361 

Xs( DPQ* , 0361 ) /POA 

Y.sl, «( ( ( , 4l+ ,  35*q4) #x)/l,  4) 

XH«<lf*H)/Ct.n.607*H) 


52 


ROO=POA*144.*XM/(53.3*TOR) 

IF<T0-68. >12/12,11 

11  FA=1,+(TQ*68# ) *1 . 85E-5 
GO  TO  13 

12  FA=1. 

13  FB=1./(1.-B4)**,5 
DO  26  1=1/3 

WDQTs359.*C*FB*FA*Y*QD*QD*SQRTF (DPD*ROQ) 

R£0=i5,28  *WDOT/(UO*QD) 

C*CO*DC*SQRTF  <10000  .  /REO) 

26  CONTINUE 
-  GO  TO  (29,40)/  1HORC 

29  T8=(Tl+T2>/2. 

TFn=<TS+TS)/2. 

TBK=CTfi*459.7>*5./9. 

TFMK=(TFM*459.7)*5,/9. 

UB=.0Q3527*T8K**1,5/(TBK+11Q.4> 

UFMs.003527*TFMK**1.5/(TFMK+110.4> 

AKFM=.001529*TFMK**.5/<l.*245.4/{TFMK*i0.**(12./TFMK))) 
RH=AC*FLL/AA 
'  G=WDOT/AC 

REFMs4,*rh*G/UFM 

REB=4.*RH*G7UB 

XCa(l.*1.915*H)/(l.+H> 

CP  =  XC*  1 241 

XN  ru  =L0GF((TS-T1)/(TS-T2>> 

U=CP*WDOT*XNTU/AA 
HS=2000 t 

XMS=SQRTF(24.*hS/(FKS*TFS>> 

YNFS  =  (TANH(XMS*FLS/12. > )/(XMS*FLS/l2. ) 
YNS=<AWS+YNFS*AFS)/AS 
HA  =  U 
DO  30  1=1/3 

XMA=SQRTF(24.*HA/<FKA*TFA) ) 

YNFA=(TANH(XMA*FLA/12« ) )/ ( XMA*FLA/12 . > 
YNO=(AWA+YNFA*AFA)/AA 

30  HA=l./< YNQ*( (l,/0)-AA*(l./(YNS*AS*HS)+TW/(AWA*WK*12. ) ) ) ) 
PR23=<UFM*CP/AKFM)**.667 

ST=HA/(G*CP) 

COL J=  ST  *PR23 

40  ALFA=  AC/AFR 

P1A=  P8*13 «  571"  PI 
P2A=  PB*13.571-  p2 
T1A=  Tl+459 ♦ 7 
T2A=  T2+459.7 
R01=.Q97318*XM*P1A/T1A 
R02=.097318*XM*P2A/T2A 

GO  TO  (41,42), IHORC 

41  R0M=.5*((TlA*R01+T2A*R02)/(TS+459.7-((T2A'TlA)/XNTU))) 

GO  TO  43 

-  42  ROAVE  =  .  097318*XM*(PlA+P2A)/(2> *TiA) 

ROM=ROAVE 

g=wdot/ac 

RH=AC*FLL/AA 

TlX=(Tl+45$.7>*5./9. 

UC=,003527*T1K**1.5/(T1K*11D,4> 

R5C=4.*RH*G/UC 
REFM=REC 

43  GO  TO  (50/45,44)/ JSQTRJN 
•44  XKD=1. 

GO  TO  60 

45  IF(R£FM-2000. >46/47/47 

A 


53 


46  T RKD-1 • 43 
GO  TO  48 

47  FM  3.049*REFM**<-,2) 

CRKD31, 09068*4, *FM*.05884*SQRTF( 4. *FM)*1. 

TRKD3i.*1.29*(CRKD-l. ) 

48  XKD=TRKD 
GO  TO  60 

50  IF(REFM  *2000 • >51,53,53 

51  SQKD  3  1.39 

52  GO  TO  56 

53  FM  3.G49*REFM**<-.2) 

54  CRKDsl. 09068*4. ^FM*,65884*SORTF< 4, *FM)+1. 

55  SQKD  3  1,*1.17*(CRKD-1. ) 

56  XXD3SQKD 

60  CC  3 .611* . 045*ALFA  * .  344* < ALFA**5 . 7 ) 

YKE3l.-2.*XKD*AlFA*ALFA**2. 

ZKC3(1.*2.*CC*(CC*CC)*( (2.*XKD)-1. ) )/(CC*CC) 
A3(ZKC-(1.*ALFA*ALFA))/R01 

83(YKE*1.*ALFA*ALFA*(1.*(4.*FD*DL/DD>  >)/R02 
FC3  (RH*ROM/FLL)*( (4,3255E*9*DPC/(G*G) )-A*8) 

PiAA  3  P1A* , 0361 
DPCAsDPC* . 0361 
GO  TO  (70#104) / IHORC 
70  GO  TO  (75/76) /  IBS 

75  COS3 .6090 
DCS3. 0410 
DSO3 .89 
BS3 . 71 
GO  TO  77 

76  COS3 .6085 
DCS3. 0250 
DSO3 . 7n 
BS3 « 56 

77  CS3 , 61 
BS43BS**4 • 

FBS31./(1.-BS4)**,5 
F AS31«  003 

DPSO= , 454*DPS  0 

PSO317,l+(TSO-220.^,376 
XS3DPSO/PSO 

Y3sl,-( ,41*.35*BS4)*(XS/1.3) 

DPSW=DPSO/.0361 
VS324.01*(TSO*218,)*.4 
USO38.25E*6+(TSO-200. )*1.45E-8 
DO  80  IS=l/3 

SDOT3359.*CS*FSB*FAS*YS*DSO*DSO*SQRTF(DPSW/VS) 
RESO3,004244*SDOT/(USO*DSO> 

CS=COS*DCS*SQRTF( 10000. /RESO) 

80  CONTINUE 

HS1831158,Q*(TSl-230.)*,5 
HS13HS18-(PS1-18.0)*.35 
HS23U55.3+(PS2A-19.  )*(  .934) 

HC3193.4*(PS2A-19, )*2.64 
WCH3WC*60 ./TC 
QAIR=WD0T*CP*(T2-T1) 

QSTM3WCH*(HSl-HC)*SDOT*(HSl-HS2) 

ERROR3 (QAIR-QSTM) *100. /QAIR 
99  PRINT  100/NR, WDOT/G/H/T1, PIAA. T2/DPCA/TS,HA/C0LU/rC/REFM/REB/ERR0R 
10QOFORMAT(I3,2X,F6.1,2X,F7.1,2X,F5.4,2X,F4.1,2X,F5,2,2X,F5.1,2X,F4,2,  .. 
.  12X/F5,1/2X/F4.1/3X/F6.5,2X/F6.5,2X/F6.1/2X/F6.1,2X/F5.1///) 

GO  TO  110  ...  ; 

104  PRINT  105,NR/WDOT/G»H/T1,P1AA/DPCA,FC/REFM 


1050FORMAT(13,2X,F6.1,2X,F7.1,2X,F5.4,2X,F4.1,2)r,F5.2/2X,F4.2»2X#F6.5# 

12X,F6.1,//) 

110  1F(NR"N0R>7/112,112  : 

112  CONTINUE 
END 
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Figure  11.  .Saturated  steam  dynamic 

vieooaity  vs.  temperature. 
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Figure  12.  Saturated  steam  specific 
volume  va.  temperature. 


Figure  14.  Superheated  steam  enthalpy 
vs*  temperature* 


. . .  -  •"*' 

67 


4 


P|Uiu*iSB _ 

»■■■«••  ■■*■»••*«!  iiimii 
■■ bbbbb "■■■■■■■■■ ■■■■*■■ 

WVlIlllllllMIIIItlllllllil 
Baa BB BBBBB BBBBB aaaBBBBaaBaaB 


aiiliiiMiim 


ssaasaBBa 


■  UBt 

six: 


■■■■■■■■■■■a 


■bbbbbbbbbi 


BBBBB BBBBBBaaaaia 


•  aaaaaaaaaiaaaaaaaaaaall 
iiaai ■iiiibiiii bbbbb ■■ 

. . . 

KMHiiiiiMiiiiiiaiiaiiM 

liiiiiiMiiiuiiBiaiijaiiii 

laaiaaaiaaaaaiiiiiaiiaiiiai 

ft" •■■«>>•«■>■■■■■■■■■ Biaaaa 
maiBiiaaaaaaiivaaaiiBtiaaa 
■  ■■■■aiaiaaaa  Biiai  al 
■aaaaaaaaBiaBaaaaaaaaaaaaal 


BBaaaBBaaaiaBaiB 


BBBBaaBBBBBBBBBBBBa 

aaaaaaaaaaaaaiaaaaa 


BaaaaaaBa _ 

■Biitiaasaa 


aaBaaaaaraa 

aaaaasBala 


MaaaaaaiiBBBiiiflB 
■iBiBaaaaiaaBBiaaaB 
BaaaiBsaaaBBaaaflBBaB 
BaBaiaiBiaiBaamaBB 


laawuaaaa 
lianafa* 
1 •••■*• »* 


’BBBBBBBBBBBBBB 


..  - - lBB_ 

■bbbbbbbbbbb 


::::::: 

■  i 


Baiiiaaa 

BBBaaaai 


Ba.aaaaaBH 


laaaaaflaaBaB 


Itaaiiaaiaan 


BaaaaaaBa 

aaaiMBH 

BBBBaaaaa 


Saturated  Liquid 
Enthalpy  vs*  Pressure 
©  KEENAN  AND  KEYES  (5) 
he  -  195.4  +  (KL9)  (2.64) 


|2  2222222222222:232224  222H1 


■■■aaaaat 


■3X1 


BBBSiaaarAaaa 
■■■■■■■aaBBaaaaaryBiBi 
Biaiatiafl aa  m  a  a a  %■  a 1 a  a  a  a | 

■•■iiaaaitaaai jmmmmmm* 
aaaaiaaaaaaaa'.taaBaaai 
aaaaaaaaaBBiriaaaaaaai 
aaiaaaaaaaaraaaBaaaBBi 
aaaaaaaaaaa <aaaBaaiaai 

BBBaaBBiaa'/BBaBBBiBaaa _ 

BBaisaaiBr aBBaaaaaaiBaaal 


■  a a bbt  aiiBBBBBiiaala 


■■•■•BBBBBBSBBBBBB 


a  b  B  a  B  B  a  a- a  ■  a 
aaiaiaaaBBB 


■BBBaaBiaa 


BBaaaaaaiBB 


flBBBBBBiaiBlia' 


laa.aaia iaBBaaBaaiia 

. . . 

3313 


BiaaBiaflaal 


■aaaar aaa 


BFjBBBBBBBBBBI 


BiaiaiaaBBBBBBaiB 


•aaaaaaaaa 


■aBaaBBaiiiiaaiaaiil 

iBBaaiaaBiBaaBBiaaBil 

BaaBaaaaaaaaaBBBBBBsl 


BBBBBBIBI BBBBIBB 


BBBIBBBIBBBI /BRla  BBII 
- BlBlaUI'.BBBBIIBa 


ViBaBBBaflBBBBBBaaaaBIBB 


iHiiaai 
liniiuiHHI 
I BBIBlIBBISlBB 
|aaaaaaaaaaaaa 

HHHBx:: 


aasaaBBaiasa 

BBBaaBBBBiaa 


aaaaaaaaaaB 


aBBBBBtaBBIBBBIB 
b  aaa  aa bbs  a  b  a  a  a  a  a 


BBBl 

i 

aril 


IBMIBBIBBB 


■aiaiiial 


133 


a a aa a  a ■ ■ aaa a a a aa a aa aaa I 
BaaaaaBBMaaBaaaaair  — 


Baiiaaaai 


BaiaflBiaiaaiaisiai 


BBBflaBBBISB 


:hH! 


•a  aaa  aaa 
aaaaaiia 


araBBBBiiBiaaiiaiaiai 


laaaBBBrjasBBBaBBaaaaB 


BIlBiaBBliBIBBI 


BBBBiaaaaiBBBaiaaiaaaa 


■aaaaaaaaBBBaaaflBaaa 


BBaaaaiBBB 


aaaaasBaaaaa 

BBBaaaBiaaaa 


BBS  ai'.BBBIB 


BaaasBBaaaBBB 


aaaiaaaaBa 

BBBBBBBBSB 

aaaiBiBlBB 


aaaiBBiaiaaa 
paaaaflaBBaaa 
Ha  a  a  aa a  a  a  a  a  a  a 

. . 

Itaiaanaapaa 

IbIbiibiubu 


Baaaaaaaia 


BBBBaaaaa 


IIBBBBBBBBBBB 


BJIBMI 

aiaaifii 


aaaaaaaaBBB 


ISBBB 

IBBBB 

iIbbb 
1  ana 
laada 

IBBBB 


3,*V>33333t 

■■■■■■riiiaBBi 
•aaaaa  /bibiiii 
BBBBI'.BiaaBBBB 
BBBSf laaBBBflaB 
■BBP4BBBS BBBBB 


aaaaaaaaa 


aaBaaBasaBaBBaaBaBBaaBaaaaal 


BBBBBBaaaBBB 


BBBaaaaaBaaBaaBBl 
■  aiiamiiiBiMi 
aaaaBBBBBaaaad 

TTTRi  b  a  a  aa  aaaaaiad 
BBBBBaBBBBBBBBaBBBaHHH 
aaaaasaaaaaaaaassaa iil 


3X3 


BBSBBBBBBBBBBrBBBBBBBB 


raaaaaaaaaaaBaaBBB 


IBBBBBBBBail 

BBBCaB8BBBBB 


b a aaaaa aaa aaaaaaaaa 


llBBBBaBBBIBBBBBBB 


BBaaaaaaaaaa 


8BBaaaaaaBaB‘,aBaBBBaaa 
aa aaaaa aaaar  aaaaa aa aaa_ 
laaaaBriBaaiBaaaBiB 


BBBBBaaiaaBB 


aaaaaaBaaaasaaaaH 


BBBBaSBBBBBB 


aaBaaaa aaaaa acv iBBaaa 


BBBBBMBBB 


BBBaaiiBiBaaaaaaBiBiaBBB, 
aaaaasaaaaaa - 1 - 


iBBBBBB 

Baaaia 


BBBl  I 

aaH) 


Baaaaai'.aaBBBBBB _ _ _ _ _ 

laaaa^aaaBiiaaiaaaiBaaBaBBBBBBiaaaBBBBSBBL 
- - 


■  BBBBBBBBIBBBB _ _ 

BsaaaaBBBaaBBBaaBBBBaaft 


aaaaaaaaaaaaa _ _ 

^ SBaaa aasaa aaa aaaaaaa aaaaa aaaaa a BBaaaa Baal 


BBS* 
a  bb  a 


!•■■■■■■■■■ BBBBBaaaBB BBBBBBBB B BBBBBB a BBSS BBBBa i  BBBB B B B BB BBBBB BBBBB BBBB B BB B B B 8 B B BB BB B B B BBBB B B B BBB B ■ BBB I 

I  33333  3;3333333333333333333333333333333333?;i:33:333:33333:3333333333333333333333333:333333333333333| 

l222222222222222!222222222522f2!!f!!!;!Sa,'..'S!2B5SB5BBB,,BBBBBBBBai,aaEBBBB,,aaaBBaBBBBa,BBB,,BBBBBBBBBB"BBl 

luilauuiaaHaiHHUHiHHiHaiaHMkiSMaiSiiSSIiiiiiiSiiaSaSHSSlMSSSMMlIXSSSSSSSSSSCSXSSSl 

SBiaaaBBBBaaaiiBaBBB 


BBBBBBBB*2JB*J  J  JJJJ  JIJJJJJBBBBBBBBBBBBBBBk-  *B B BB BB BBB B BBBB BBBPB BBBBB BBBBB BBBB B BBBB ■ 

I:::::::::::::::::::::::::: . 3x3333333 . 

laBBBaaaaaaaaaaaaaaaB.aaaBiBaBaaaaaariBaBaasiBaaiiasasBaafl] 


aaaBaaBaBaBBBaBaaaaBa BaaaBaaBaaaaaBBaaaaaaBBaBBaiiiaaiaiaaaBi 


aaaaa aaaaBBBaaaaa 


lUBiiBRaiiBBBBBiBaaaBBaiiaiai 


a  aaa 
•aaa 
Via 
aa  aa 
a  aaa 
aasa 


12222522222222222222222222222222222222222 222B2fBBBBB"BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 
I2222222222222222222222222222222-2222222222222bbbb!bbbbbbbbbbbbbbbbbbbb!,,bb5bbbbbbbb*bbbbbbb,,b,,bbbbbbi 
^  ;  :  :  . . . 


|BaaaaasB<- 
| IBBBBIBBI. 

■  •BBaaBBBaaaaaBBBBaac'.B _ _ _ _ _ _ 

(5«aaaBBaaaBaaaBBBBBr aaaaa aa aaaaa aaaaa aaa  aaaaa aaaaa aaaaaaaaa a 


aiaaBBiaBBsiBi a aaaaa ■ 


bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbsH 

aBBBBBBBBBBBflaBBBIBBBBIBBBflBaaail 

aaaaaaaaaaaaaasasaaBaBBaaaaaaaaaaaaal 


■  «*■■■■■  a  ■*«*>/' *fc«Baaa] 
siaxaHaaaaai^^aaasBaal 
aaaaa  aa  aaa  ai'>.ai  8  aaaaa 

^^Hmaaaanaar  aaaaa  aaialHHHHHiHHVHHVHHHHVHHOTI 

. . . 

•■•■■•■■afliaa.aaar ja aaaaa aaaaa aaaaa aaaaa aaaaaaaaaaiH 
BBBBB  a  B  aaa  I  a  aaa  a  a  .as  is  Baa  BBBBaaaaa  a  aaaai  iia  aaaaa  iraH 

a  b  b  b  a  a  a  a  aa  a  a  a  a  a  a  1  a  1  a  a  a  1  a  a  a  aeiaiaaiiTaiiaT] 

■  ■•■■Baa  mb  aaaaa  iiaiaaaaajraaaaasaaaftaaaaal 
a  a  a  a  a  a  nAl  a  ■  a  a  a  ■  a  a  a  a  a  a  a  a  a  a  a  a  aa  a  a  a  aft  a  a  a  a  a 

|BBBBBamfBBBBB«aBBBBBBBBBBBBlaBBBIBBBBBBBM 

■HMaBaaiaaiBaBaaBiaaBsaBaaiBBBflaaaaaBaaaBBaBlBBBBBiBiiBflBaaBBBaBaBifliaBBa 
BiaBiiiiaaiaiaBaaiaaiiaiiBaaii aaa aaa a  a aaaBaaa aaa aaaaaaaaaaa bib bib aaa as 
■■■■■■ ■■■■■•aaa a a aaa aaa as aaaaa a aaa a aaa a a aaa aaaaaaaaaaa a a aaaaa a aaa aaaaa 


■■■■■■■ ■■■!■■■■■■■■■■ aaVaas aaa ■■■■■■■■■■■■■!■■ ■ I 

■  aaBBaflBBBBaailBBBBBBBBBBBBBBBBBB jBBBBIBBBB SB B  B  B | 

- BSBBBaaaaaaBaaBSBBBBBaaBaaaaaaaaaaBl 

__  (BiBliBaBBiiBBiiaBBaaaBiaBBaiiiiiaBl 

r>i~ - ^3  8  B  ■  B  ■  8  BB  BB  BBBBB  BBBB  B  BB  B  B8  BB  8  B  ap  B  B  B  B  B  B  ■ 

iBBBBaBBBBBBBBBBBBBBBBBBBBBBBBBBBBBaBBBBBBBBBBBB  BBS  B B  BBB  BB  BaB ■ B  BBS B  B | 


laBaBBBBBnUIBBB'laiBlllaiBBS 
!■■■•■  •■•■■■■■Bar  XaaaBBBBBaBMJ 
(• . . 

|  a  a  b  b  ■  b  a  v  a  a  a  a  a  a- ,  a  a  a  ■  a  b  a  b  a  IHHHHHHVHIVHHHHVHHVHHHHHVMHVHHHiWHHHHnftnViWVHWlHHHWVnVVW 

1 222222  2222222  j2  22222  222252222222552S2222iMMlMSS£SSSIlHllMHHHHHHHHHHHHHHHHSHaHHSHHMH£HMS£MS£2H 

ft ■•■■■•  ■•••■*  ^aa ••■•■■■■■■■  aaa aBMMMMgJ 


IBBBBB;B' §!■■•■■;■■■■■"»■■■■■■■■■*■■■•■■■ ■■■■■•■■■BBBBBBBBBBB ■■■■■■■■BBBBBBBBBBBBBBBBBBBBBBBBBBaaMaBl 
••■■■■*■■■■■■■•■■■■■■■■■■•■■■■■«■■■■ ■■■••■■■■•■■■BBBBSBBBBBBBBBBBBBBBBBBBBBBBBBBBBBaa BBBBB BaBBBBBBaBl 
■§■•■•*■■■■■■■■■■■■■■■■■■■■■■■*■■■■■■■■■ ■■BBBBBBaBaBBaBBaBBBaBBBaBBBBBBBaBBaBBBBaBBBaBBBBBBBaBaBBBasl 

2222?^2222222222!2;22222bbbbbbbbbbbbbbbbb-bbbbbbbbbbbbbbbbbbbb^bbbbbbbbbbbbbbbbbbbbbbbbbbbb»"b"»»»b"I 

BBaar  BBaaa  ^a  bbbibbbb  a  aaa  bbbbbb  bbbbb  Baaaa  saa  a b  bb  a  aa  bb ■■ bbbb  aa  sa  aaa  aaaaa  as ■■ a a  a  aaa  aaa  aaaa  aaa  aaaaa  a  a  aaa  | 

!!2f.2!!!22!2!!!!!!!!!!!aaaailBailliaa*BBtBaBBBBaaaaBaaaa**aa*a**IBa*aaitaBaaaaaaaaaaiaaBaiaiBa*"""BBa""Bi 

. . . . a . . 

■■^•■■■■■g ■■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■aBBaBBBBBaBBrBBBiBBBaaaBBaBaaBBBSaBBBBBBBBBBBBaBBBBaa 
F<2!!!!!!!!!!!!!iaBI!aaaaaaaBaaB"BmBBaBaBiBaaaBa8aiaBaaBaBBaaaaiaaa|aaaBBaaaaiBaaBBBaBBiiBBBBBaaBB 
'^•■* ■•■■■■■■■■•■■■■■■■■• ■ ■■■■■■ ■ ■■■■■■•■■■«■ ■■ ■■■■■■■ a ■••■ ■ aaa a aaa aaa a a ■■■ a aa as saa ■■■■■■■ saa BBBBaai 

■■■■■■■■■■■■■■BBBBBaaaaBBBBBBBaSBBBBBaBBaaaaaBBBaBaBaBBBBBBBBaBBBB 

2r22222!f222!!BBB22B!!BBBBBBBBBBBBB,r*B,,BBBBBBBBBBBBBBBBBBBf  ■  ■■■■■■■  BBBBBBBBBiBBB  ■■•■■■•■■■■■■■■■*■■■ 

2222B2BBBaB2BBBBaaaa  ■■*•■*■■■■*■  ■■■■■■■■  BaaaaaaaaaaiaaaaaaM  aaaaaaaaaaaaaaaaaaaa  ■aaaaaaaaaaaaaaaaMa 
BBBB  ■■■■■■  ■  ■■■•■■■■■■■■■■■••■■  bmbb  aa  aaa  ■■■■■■  ■  aaa  bbs  BaaBaaBBt  ■  ■■■■■«■■■■■■■■■  as  BBBaBaaBaBaaBBBaaBBB 
222222222222222222222BBBaBBBBBBBa9£BBBBB ■■■■■•■•■■■■■■■■■■■■■■■■■•■■■■••■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■ 
2222g22222BBBBB,BBB-,*BBBBBBBBB ■■■■*• ■*■■■■■■§■■■■■■■■■■■■■■•■■■•■■■■■•*■■■*■■■■■•■*■■•■■■■■■■■■■■■•* 
■■■■■■■■■• ■ ■■■■■■■■• mb . . ibbbb a si BBBaa as iBasaaBBa a bbbbb ■■■■■ iibb ■■■■■■ bbbbb bbbb fl 

■■aaaBBaaBaaaaBBaaaaaaBaBBaBaaaaBBBBBBBBBBaaaa.aaaaaaaaaaaBBaBBaaaaaaaaaaaaaaaaaBaaaaaBaaaBaaaBaaaaas 

■■■■•■■■■•■■■••"■•••BBBBBBBBBBBBBBB  B  B  B  B ■ SB ■ SB  BBS  B  B  8  BBBBBBBBB  BBBBBBBaaaS  BB  B  BBBBBB  BBBBBBBBBB  B  B  B  BS  B  B  BBS 

*■■•■■■■■■■■■■■■■■■■■■■■■■ ■■•■■■■■■■■aBBBaaaa bbbbb BaBaaaaBaaBaaaaBBBaasaaaaBBaaaaaaaaBBaaBBaBB-aBaaaa 
222222222222BBB2BBBBBBBBBBBBfcBBBBBBBBBBBBBBBBBBBBB •■■■■■■■■■*■■■■■■■■■■■■■■■•■■■■■■■•■■■■■■■■■■■■•■■ 
2222222BBaBaBaaBBaBBBBBBBBBBBBBB,BJBB  ■•■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■•■■■■■■■■■■■  ■■■•■■■■■■■■•■■■•■■■ 
■■■■■■•■■•■■•^••■■■a ■••■■•■■••■■■■■■■■■■ ■■■■■■■■■■■■■■■■BW»BBBaBWfB»BBBBBBBBBaBBBB«BBBaBBBBaaB BBBBB 


